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ABSTRACT 

We present mid and far IR ISO data of the CfA Seyfert galaxy sample. These data allow a 
detailed study of the far IR Spectral Energy Distribution (SED) of these galaxies. A Bayesian 
inversion method has been used to invert the SED of these sources yielding two fundamental 
results, namely, that the mid and far IR SED of Seyfert galaxies can be explained solely through 
thermal reradiation of high energy photons by dust, and that this thermal emission is made 
up of two or three different independent components, a warm, a cold and a very cold dust 
component. These thermal components have been readily explained as produced respectively by 
warm dust heated by either the active nucleus or by circumnuclear starburts, cold dust heated 
by star forming region in the galaxy disk, and very cold dust heated by the general interstellar 
radiation field. Comparisons between the parameters obtained from the analysis of the IR SEDs 
(fluxes, temperatures, luminosities) have been made. Our results suggest that the emission in 
the mid IR is anisotropic and the differences found between Seyfert 1 and Seyfert 2s can be 
explained with thin molecular tori models. 



- 3 - 



1. Introduction 

The origin of the IR emission from Seyfert galaxies has been a matter of discussion since the early work 
of Rieke & Low (1972); Stein (1975); Rieke (1978); Neugebauer (1978), and others. These authors found an 
emission excess between 3 and 5 /im in Seyfert galaxies, that originated a strong controversy between those 
who supported that the emission was non thermal and those who defended that the excess was emission 
from dust heated by the nucleus. The use of efficient instruments to observe at 10 /im, and especially 
the pioneering work at the far IR carried out from the Kuiper Airborn Observatory (KAO) showed the 
importance of the mid and far IR emission to quantify the bolometric luminosity of Seyfert galaxies (Telesco 
& Harper 1980; Smith et al. 1983). However, it was the IRAS satellite that provided an extensive set of IR 
data for a large number of galaxies, from which it was shown that Seyfert galaxies are indeed strong far IR 
emitters (Rodriguez Espinosa, Rudy & Jones 1987; Edclson, Malkan & Rieke 1987; Spinoglio ct al. 1995). 

The IRAS satellite was key to the understanding of the importance of the IR emission to the total 
luminosity of active galaxies. However the IRAS data alone are not sufficient to clarify the nature of 
the IR emission, as there are measurements only at a limited number of wavebands preventing a good 
definition of the shape of the Spectral Energy Distribution (SED) at mid and far IR wavelength (see, e.g., 
Telesco 1988; 1990; Bregman 1990). A proper characterization of the mid and far SED is essential to 
understand the emission mechanisms that produce the high output of Seyfert galaxies in the IR. Recent 
studies have discussed the origin of the IR emission suggesting that it is of thermal origin. For example 
Giuricin, Mardirossian & Mezzetti (1995) have studied a complete sample of Seyfert galaxies at 10 /zm and 
propose that the emission is due to thermal reemission by dust. Bonatto & Pastoriza (1997), based on color 
studies of IRAS data from diverse Seyfert samples, find that the colors obtained can be explained with a 
combination of dust heated by the nucleus plus cold dust in the host galaxy. Siebenmorgen et al. (1997) 
show that the SED at IR and milimetric wavelengths can be modeled assuming that the dust heated by 
a central source dominates the luminosity output of these objects. Maiolino et al. (1998) confirmed this 
last result finding, in their high resolution IR images of the Circinus galaxy, an unresolved source, with 
size < 1 pc, that is reprocessing the nuclear output via dust reradiation. Rigopoulou et al. (1997) have 
observed a sample of AGN in the CO mm line suggesting that the far IR emission of Seyfert galaxies is 
thermal, based on three different evidences: the correlation found between the far IR and the CO emission, 
the dependence of the far IR emission to hydrogen molecular mass ratio with dust temperature, and the 
similarity of the profile shape of the CO and HI lines. 
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Another important issue is the understanding of the differences between the two Seyfert types. 
According to the unified models, Seyfert 2 nuclei are intrinsically similar to Seyfert 1 nuclei, the differences 
observed being due solely to geometrical effects. In Seyfert 2, neither the broad line region nor the optical, 
UV and soft X ray continuum can be observed directly because the central region is obscured by intervening 
material in the line of sight. Some authors argue that this material forms a sort of disc or torus of molecular 
material. This disc or torus is thought to be reponsible for the collimation of radiation and the observed 
anisotropics, i.e., biconic structures in emission line images (Simpson et al. 1997; Wilson et al. 1993) or 
highly collimated jets. It is expected that at sufficiently long wavelengths the optical depth of the torus 
would decrease and the differences between the two Seyfert types should disappear. Several tests have 
been made to ascertain the presence of these molecular tori (Hcckman 1995; Pier & Krolik 1993). These 
obscuring tori have been theoretically modeled by Pier & Krolik (1992) and Granato & Danese (1994) 
among others, predicting that the mid IR optical depth is still considerable, thus it should be expected that 
Seyfert Is are more luminous than Seyfert 2s in the mid IR. 

In this work, we make an attempt to understanding the origin of the mid and far IR emission from 
Seyfert galaxies through the study of their spectral energy distributions (SED). We present ISO data of 
the entire CfA Seyfert sample (Huchra & Burg 1992), consisting of 25 Seyfert 1 and 22 Seyfert 2, plus a 
LINER. Section 2 describes the observations. Section 3 presents the separation of the Spectral Energy 
Distributions in thermal components by means of the Inverse Planckian Transform. Section 4 describes the 
thermal emission components obtained from the inversion and in Section 5 we perform a statistical analysis 
of the parameters obtained and discuss the differences between the two Seyfert types. 

As the CfA Seyfert sample is a complete sample of Seyfert galaxies we expect that the results obtained 
here are statistically significant for all Seyfert galaxies and certainly suggestive for other classes of AGNs. 

2. The data 

Observations of the CfA Seyfert sample have been carried out with the Infrared Space Observatory 
(ISO; Kessler et al. 1996) through filters at 16, 25, 60, 90 120, 135, 180 and 200^m. NGC 1068 was also 
observed at 4.9, 7.3 and 11.5 /im. The filter set was chosen to achieve good coverage in wavelength while 
producing a good sampling of the SEDs. The ISOPHOT instrument (Lemke et al. 1996) was used in 
the PHT-P and PHT-C configurations, with the PI, P2, C100 and C200 detectors. PI and P2 are single 
element photodiodes. C100 is an array of 3x3 pixels, each one projecting onto 45 arcsec in the sky, while 
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C200 is a 2x2 pixel array, each pixel projecting onto 89.4 arcsec in the sky. Integration times were calculated 
from the signal/noise ratio estimates produced by the ISOPHOT simulator based on interpolations and 
extrapolations of the IRAS data. The PHT-P observations were done in chopping mode through a 120 or 
180 arcsec aperture depending on the size of the objects. The PHT-C observations were done in staring 
mode. In this case, to set the background level, an empty area adjacent to the source was observed prior to 
the actual source and with identical instrumental settings. For NGC 1068 we have done maps with C100, 
at 60 and 105 /im, moving by half a pixel (23") between two contiguous map positions. Four additional 
objects (NGC 3079, NGC 3227, NGC 4051 and NGC 4151) have been mapped at 90 /jm with the C100 
array. Details of the mapping observations of these four galaxies are given in Perez Garcia, Rodriguez 
Espinosa & Fuensalida (2000). 

The data reduction was done with the help of the PHT-Interactive Analysis (PIA) V.7.0 tool. PIA 
first deglitchs the data eliminating the cosmic rays and others spurious effects, then corrects for non-linear 
effects taking into account both the dynamic range of the detectors and the flux level of the observed 
sources. To determine the signal PIA linearizes the integration ramps, and corrects for detector drifts. 
The background substraction is achieved by repeating the PIA reduction process for the background 
measurements files, which are then substracted from the object files. It is also necessary to correct for the 
PSF fraction that is not seen by the detector (see Lemke et al. 1996). This factor varies with wavelength, 
hence the actual correction factor for each filter is built-in in a PIA table. The flux calibration is performed 
using the standard calibration of PIA 7.0. For the PHT-P measurements the process is similar except 
that the background substraction is performed automatically from the chopped measurements and the flux 
calibration is done through the use of the internal FCS (Fine Calibration Source; Lemke et al. 1996) since 
it provides a better match with the actual intensity of the sources. Note that the P detector response is 
very dependent on the intensity of the source being observed. 

The photometric uncertainties are small varying between 1% for the objects with higher signal/noise 
ratio and 15% for the weakest objects. The final uncertainties are however dominated by the uncertainties 
in the flux calibration. The calibration factors have kept improving as the detector and instrument 
characteristics have become better known, although there are still some residual effects (non linearity of 
detectors) that are yet not adequately modelled. As of this writting, the uncertainty of the calibration is 
equal or better than 30%. Therefore, throughout this paper, we adopt a conservative uncertainty of 30%. 

The fluxes for the observed objects are given in Table | and |. Fluxes are in Janskys. Some objects 
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could not be observed in some or all of the niters for various reasons. These are also given in Table 1. 

Figure [I] shows the Spectral Energy Distributions (SEDs) of the observed objects. The ISO data have 
been plotted together with IRAS data (Edelson, Malkan & Rieke 1987) at 12, 25, 60 and 100/mi showing 
that the agreement of the ISO and IRAS data is very good in most cases. 

EDITOR: PLACE TABLE | HERE. 
EDITOR: PLACE TABLE | HERE. 

3. The inversion of the SEDs 

It can be seen in fig. 1 that the SEDs describe a well defined energy range with a steady increase from 
the near IR on to a maximum between 90 and 135/im and the start of a decline toward longer wavelengths. 
Several results like those mentioned in the introduction point to the mid and far-IR emission of Seyfert 
galaxies being of thermal origin. Furthermore, we have shown that the SED of a few Seyfert galaxies can be 
explained as emission produced by two emissivity weighted blackbodies (Rodriguez Espinosa et al. 1996). 

EDITOR: PLACE FIGURE @ HERE. 

However, rather than proceeding with a plain fit of a "ad hoc"number of black bodies, we have 
preferred to use an inversion method to analize the SEDs of the galaxies in the sample. This has the 
advantage that no assumptions have to be made as to the number or location of the sources responsible 
for the observed spectrum. In particular, we have used an Inverse Planckian Transform algorithm, that 
employs an emissivity (e oc A" 1 ' 5 ) weighted Planck function kernel to switch from frequency space to the 
temperature domain, hence revealing the temperature distribution of the sources that originate the observed 
SEDs. The method applies Bayes theorem of conditional probability and the Richardson-Lucy iteration 
algorithm which converges quickly to a optimum result. To increase the number of data points used in the 
inversion algorithm, whenever available, we have added the four IRAS band fluxes at 12, 25, 60 and 100 /im 
as given by Edelson, Malkan & Rieke (1987). Furthermore, to avoid boundary convergence problems in 
the inversion algorithm we have used 10 /jm ground-based data from several authors (Rieke 1978; Edelson, 
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Malkan & Rieke 1987; Maiolino et al 1995), and 1.3 mm upper limit data from Edelson, Malkan & Rieke 
(1987). These additional data have been used solely for the purpose of constraining the inversion algorithm 
at the borders of the wavelength range of interest. Details of this method are given in Salas (1992) and in 
Perez Garcia, Rodriguez Espinosa, & Santolaya Rey (1998) and are not repeated here. 

Figure [2] shows the results obtained after application of the Inverse Planckian Transform to the SEDs 
of the CfA Seyfert galaxies observed. For each object, the upper pannel in Figure || plots the ISO data 
(triangles) and the IRAS data (filled squares) and the best fit to the mid and far IR data in heavy black, 
with the different spectral components contributing to the fit printed in dashed lines; the bottom pannel 
shows the temperature spectrum that produces these components. 

EDITOR: PLACE FIGURE | HERE. 

It is important to note that the temperature spectra obtained in this way are continuous. Further, the 
spectral temperatures of the dust grains group themselves in discrete and well defined features that we have 
called components. These components can be described by their peak temperature. Figure || shows that 
34 out of the 40 galaxies for which we have obtained good temperature spectra show three temperature 
components: a warm component with central temperature T~150 K, a cold component with T^40-70 K, 
and a colder component with T~15-25 K. The remaining 6 objects show only two components, although 
in NGC 1068 when the short wavelength data are added a third component is recovered. Table || shows 
the peak temperatures of these thermal components for each object. This table also gives the spectral and 
morphological type for each object. These results confirm the analysis made before for a few objects of the 
same sample (Perez Garcia, Rodriguez Espinosa, & Santolaya Rey 1998). 

The fits obtained from the inversion are, in general, very good. In some cases, weak objects or objects 
with low signal to noise ratios, the uncertainty in the inversion is higher. For some objects, we have 
rejected some measurements (the flux at 120/J,m of NGC 4235, and the fluxes at 90^m of Mrk 270, Mrk 461 
and 1614+35) because these data points do not follow the shape of the SEDs defined by the rest of the 
measurements. In all of these cases the rejected measurement is very noisy. Four other objects (Mrk 590, 
Mrk 573 and Mrk 841) have not been measured at long wavelengths (>90^m), because they are below the 
sensitivity limit for detection with ISO. 

Nevertheless, the method still finds three components for three out of the except for NGC 1068, that 
will be analized in a separate section, four objects. For Mrk 590, the inversion does not converge. 
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In four of the objects, the warm component peak temperature obtained from the inversion is not high 
enough to fit the data. These four objects are Mrk 335, NGC 5273, Mrk 1243 and IZwl. Both in Mrk 335 
and in NGC 5273 it can be due to the low signal to noise ratio of the data. NGC 5273 shows a discrepancy 
between the IRAS value at 12 /im and the ISO value at 16 /jm (again very likely due to the low signal to 
noise ratio of the ISO data point) that prevents a good fit to the data. Mrk 1243 has only three data points 
in the hot part of the spectrum. IZwl is a special case, since the short wavelength data are of good signal to 
noise ratio, however it is not detected at 200 /j,m and the flux at 180 /im is rather low. The fall of the signal 
at longer wavelengths is fast enough to explain that the object is not detected at 200 fj.ni. Figure || shows 
that the SED of IZwl if made of just two components does not account for all of the observed emission. If 
we substract this initial two-component fit from the actual SED, the residual seems to correspond to yet 
one third hotter thermal component. The inversion method does not succeed at identifying completely this 
third component, due to the temperature of this third component being too high, therefore only a small 
fraction of the energy within this third component falls within the range of wavelength studied. The shape 
of the residual is however similar to the other two components, a clear hint of the cxistance of yet another 
thermal component with a peak temperature of ^300 K. This temperature is the highest temperature found 
among all of the SEDs. The cold and very cold components also have peak temperatures higher than the 
corresponding components of the SEDs of the rest of the CfA sample. IZwl is one of the more distant and 
brighter Seyfert galaxies known. Its different spectroscopic characteristics of line width and line ratios were 
pointed out by Osterbrock and Pogge (1985) and Goodrich (1989). Recently, several authors have pointed 
out the different nature of IZwl. For example, the optical and X-ray spectra are steeper in IZwl than 
in normal Seyfert 1 and 2. Indeed IZwl shows fast and large amplitude variability in its X-ray emission 
(see Halpern & Moran 1998 and references there in), and is possibly better associated with the broad line 
quasars (Boroson & Meyers 1992; Lawrence et al. 1997). 

4. The dust components 

The inversion method assigns to each dust component a range of temperature, i.e., the dust grains 
have temperatures ranging between a minimum and a maximum temperature. This is indeed an expected 
and physically meaningful result, as the dust closer to the source or sources of radiation will be warmer 
than the dust farther away. The existence of temperature components indicates that there are well defined 
sources with well defined temperature profiles which must be the result of well defined physical scenarios 
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existing in these objects. From the temperature spectra, we can calculate the flux that each dust component 
contributes to the total mid and far IR emission. For each component, the flux is obtained by integration 
of the temperature spectrum over the relevant range of temperatures: 

r-Trnaxi 

Fi= ^{T)dT i=l,2,3 (1) 

J Tmirii 

where T m i ni and T ma:Ei are the extreme values of the temperature range that defines each component, 
and \1/(T) is the temperature distribution obtained from the inversion process. Table ^| shows the fluxes 
obtained in this way for each component and the total flux for each object. The comparison between the 
far IR IRAS fluxes (FIR), as calculated from the values at 60 and 100 /im (Lonsdale et al. 1984) and the 
sum of the cold and very cold components calculated here is in good agreement. The luminosities for each 
of the components are shown in Table ||. We have used Hq — 75 Km s _1 Mpc -1 throughout this paper. 

The key point now is wheather these emission components can be physically explained within a sensible 
scenario. In what follows we review the pieces of evidence that we have for explaining each of components: 

a) In a previous paper Rodriguez Espinosa & Perez Garcia (1997) were able to use optical R band 
images of a subset of low redshift Seyferts of the present sample to separate the fluxes from the central 
region from those of the galaxy disks. We found a very good correlation between the ratio of the extended 
to the compact R band fluxes and the ratio of the cold plus very cold component to the warm emission 
component. The conclusion was that the warm emission component is related with the central regions of 
the galaxies while the cold and the very cold emission had to originate in the disk of these galaxies. 

b) Furthermore, a correlation has been found between the flux produced by the warm component and 
the flux in the high ionization coronal lines fluxes like [OIV]A25.9 and [NeV]A14.3/im (Prieto, Perez Garcia 
& Rodriguez Espinosa 2000), indicating that the warm component must be heated by the nucleus of these 
galaxies. 

c) In another recent paper, Perez Garcia et al. (2000) have shown that based on 90/im ISO maps of 
four nearby Seyferts, the 90/xm emission is physically extended up to radii similar or larger than those 
seen in optical images of these same galaxies. Furthermore, the extension of this 90/xm emission has been 
characterized and its extension, scale length and surface brightness profiles are typical of normal galaxy 
disks. 

Based on the above, a scenario arises in which the warm component is associated with dust heated by 
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radiation coming from the nuclear or circumnuclear regions of these galaxies, while the cold and very cold 
dust must be heated by process occurring in the galaxy disk. Danese et al. (1992) also conclude that the 
mid IR emission (10-25/xm) is dominated by the nucleus or the circumnuclear region. Further support to 
this scenario comes from the following: 

• Warm dust. Its characteristic peak temperature is in the range 120-170 K, a range of temperatures 
warmer tham is normal of dust in typical starforming regions. The nuclear origin of the radiation 
responsible for the heating of the warm dust was already indicated by Rudy (1984) who found 

a correlation between the [OIIIJA5007 emission line flux and the 10 /xm emission for a sample of 
quasars, Seyfert galaxies and radiogalaxies. This author suggested that the dust responsible for the 
10 /zm emission is mixed with the ionized gas of the Narrow Line Region (NLR) that produces the 
[OIII]A5007 emission line. More recently, Giuricin, Mardiossian and Mezzetti (1995) support the idea 
of the nuclear heating of this warm dust, based on 10 jim small aperture observations of a sample of 
100 galaxies. They found that the 10 /zm emission correlates very well with the 25 /an IRAS emission, 
while the correlation is poor with the 60 /xm emission. From a different perspective Heckmann et al. 
(1997) have found in the Seyfert 2 galaxy Mrk 477 a strong starburst with very warm dust in an scale 
of a few hundred parsecs, that confirms the idea that the warm dust can also be heated by nuclear 
and circumnuclear starformation regions. 

• Cold dust. Peak temperatures for this component range between 40 and 70 K, a range of temperatures 
that is typical of dust in regions of starformation. Note that cold dust is present in mostly all classes 
of galaxies, including normal and starburst galaxies (Knapp et al. 1996; Chini, Kruger & Kreysa 
1992; Klaas et al 1997; Walterbos & Schwering 1987). In all these galaxy types, the heating of the 
dust is produced by OB stars in star formation regions in the galaxy disks. The higher temperatures 
corresponding to higher recent starformation rates (Young et al. 1989) 

• Very cold dust. Peak temperatures for this component range from around 15 to 25 K, temperatures 
that are normal of dust heated by the general interstellar radiation field. This very cold dust is 
tipically composed of big grains in thermal equillibrium with the interstellar radiation field and 
has been observed in normal spiral galaxies (see, e.g., Walterbos & Greenewalt 1996; Walterbos & 
Schwering 1987; Cox, Kruger & Metzer 1986). For example, Cox, Kruger & Mezger (1986) predict 
that the dust of the Galaxy is at a temperature of around 15-20 K, and Walterbos & Schwering show 
that in the disk of M31 there is very cold dust at 21 K. 



4.1. NGC 1068 



NGC 1068 has been observed at a wider spectral range than the rest of the sample, 11 filters between 
4 and 200 (im have been employed. NGC 1068 is indeed one of the best studied Seyfert 2 galaxies in all 
spectral ranges. Morphologically it is clasified as SAb, and was the first galaxy in which broad emission 
lines in polarized light were found (Miller & Antonucci 1983; Antonucci and Miller 1985). In the IR range, 
Tresh-Frenberg et al. (1987) have found that the emission between 5 and 20 /Ltm presents an elongated 
structure in the nucleus, that they explain as thermal emission by dust heated by both the active nucleus 
and an active starforming region near (l"-2") the nucleus. The comparison between 10/im images and HST 
images (Cameron et al. 1993) indicates a correlation between the hot dust emission and the Narrow Line 
Region. These autors developed a model that replaces the molecular torus by giant molecular clouds of 
gas and dust that attenuate the broad lines. Braatz et al. (1993) found a correlation between a 12.4 /im 
image and both optical continuum and [OIII] images; they also observed that the IR emission is aligned 
perpendicular to the torus plane and infered that the extended IR emission is thermal radiation by dust 
in molecular clouds heated by collimated nuclear emission, without ruling out a possible contribution by 
dust in the molecular torus. Telesco & Harper (1980) have studied the far IR emission of NGC 1068 with 
KAO (Kuipern Airborn Observatory) observations between 30 and 300 fim and have found two thermal 
components with temperatures of ~36 K and ^115 K. 

We have applied the inversion method explained before in the same spectral range that we have 
used for the rest of the objects in the sample, i.e., between 12 and 200 /im. Figure || shows the results 
obtained, i.e., two thermal components, a cold one with peak temperature of T^36 K and a warm one with 
T^115 K in the peak. It remains however to extend the fit to the shorter wavelength data which we have 
for NGC 1068 alone. When we apply the inversion method to the full range observed (4-200 ^m) the SED 
of NGC 1068 appears now as the sum of three thermal components (see Figure ||). Two of them can be 
identified with the two thermal components previously found before without considering the 4-12 /im range. 
The peak temperatures of the three components are now 29, 110 and 278 K. The differences between the 
flux contained in each of the two components common to both inversions is less than 7%. 

EDITOR: PLACE FIGURE | HERE. 

The three components found in NGC 1068 are explained as in the rest of the objects of the CfA sample. 
The cold component is produced by dust in the host galaxy, the warm component corresponds to dust 
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heated by the active nucleus and/or circumnuclear starforming regions. The dust at higher temperature is 
also heated by the active nucleus, and should either be located in the inner part of the torus or be heated 
by active starforming regions placed in a radius <100 pc from the active nucleus (Gonzalez Delgado et al. 
1998). The flux enclosed under the component with T-278 K is F = 9.17 1CT 9 erg cm" 2 s -1 , and the 
corresponding luminosity L = 2.40 1 45 erg s -1 , while the total IR luminosity of NGC 1068 is 8.60 10 45 erg 
s -1 . Hence this component represents 28% of the total IR luminosity, which agrees with the fraction of the 
total UV emission that is produced by circumnuclear starforming regions in a sample of Seyfert 2 galaxies 
(Gonzalez Delgado et al. 1998). 

5. Energy balance: The cold emission component 

One possible test that can confirm the suggestion made before on the relation between the cold emission 
component and the presence of starforming regions in the galaxy discs consists of testing the ratio of Ha 
emission to far IR output. For instance in NGC 3079, the observed Ha luminosity is 7.93 10 40 erg sec -1 
(Armus, Hcckman & Miley 1990), a figure that includes both the nuclear and extended Ha luminosity. 
Moreover, with a ratio of extended to nuclear Ha emission of 0.77 (Armus, Heckman & Miley 1990), the 
extended Ha emission amounts to 3.49 1 40 erg sec -1 , which shows the importance of the Ha flux produced 
in star forming regions outside the central region in this object. If we take as an average extinction the 
value measured by Hawarden et al. (1995) Ay =7. 5 mag, the extinction at H a is K Ha =&.l mag, and the 
extinction corrected extended Ha luminosity is 0.93 1 43 erg sec" 1 . This is to be compared with the IR 
output due to dust in starforming regions that competes with the gas for high energy photons from massive 
stars. In NGC 3079 this is what we have called the cold emission (Perez Garcia, Rodriguez Espinosa & 
Santolaya Rey 1998), which amounts to 3.52 10 44 erg sec -1 , or a ratio of Ha to cold FIR emission of 0.03. 
This value is in agreement with the expected value in regions where massive stars are forming (Devereux 
& Young 1990). Hence for NGC 3079 there is good agreement between the Ha output from disk HII 
regions and the cold IR emission, which confirms the results in Perez Garcia et al. (1998) and in this work 
indicating that the cold component emission in Seyfert galaxies originates in starforming regions wjthin 
their disk. 

In the case of NGC 4051 we consider the ratio between Ha and far IR luminosity in two limiting cases. 
First, we consider the total integrated Ha luminosity, 4.34 10 41 erg sec -1 (Romanishin 1990), which after 
correction for extinction (Ay= 0.24, Ho et al. 1997) becomes 5.22 10 41 erg sec -1 . Hence, the ratio of total 
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Ha to cold far IR emission (0.25 10 44 erg sec -1 ) is 0.02. This is an upper limit as the total Ha emission 
includes the nuclear emission which to a certain extent can be due to non-thermal processes associated with 
the AGN. On the other hand, if we consider the Ha emission coming exclusively from the star-forming 
regions in the disk of NGC 4051, which according to Gonzalez Delgado et al. (1997) is 26% of the total 
emission we get a ratio of extended Ha to cold far IR emission of 0.006. This is a lower limit as some of the 
Ha flux that we have assumed nuclear is very likely due to extended circumnuclear emission that will also 
contribute to the cold far IR emission. The actual value for the ratio of extended Ha flux to cold far IR 
emission ranges therefore between 0.006 and 0.02, a value that is well within the range given by Devereux 
& Young (1990) as typical of normal HII regions in the disk of spiral galaxies. 

In the case of NGC 5033 the integrated Ha flux is 3.98 10 41 erg sec" 1 (Kcnnicutt 1983), which after 
correcting for extinction (Ay= 1.48, Ho et al. 1997) becomes 1.20 1 42 erg sec -1 . The Ha emission to 
cold far IR emission ratio is 0.02, which is again to be taken as an upper limit given that the total Ha 
flux includes the nuclear emission which we can not separate because we do not have the necessary data, 
However this upper limit is similar to that obtained for NGC 4051. 

As for the rest of the sample a similar exercise with Ha fluxes obtained from the literature produce 
similar results. In very few instances have we been able to separate the extended from the nuclear Ha 
emission. For these cases the ratio of extended Ha emission to cold far IR emission ranges between 1. 10~ 5 
for NGC 5929 and 0.02 for Mrk 270. In most cases this separation between extended and nuclear emission 
has not been possible and hence ratios between the total Ha emission and the cold far IR output have 
been obtained. These range between 1.5 10~ 4 (Mrk 766) and 2. 10~ 3 (NGC 3516). In all cases it can be 
seen that the ratios are within those expected if massive stars are responsible for the ionizing photons as 
well as for the heating of the dust. We therefore conclude that the cold emission component that appeared 
naturally from the inversion process is related to radiation from dust in starforming regions in the galaxy 
discs. 

6. Seyfert type differences 

The current understanding of Seyfert galaxies assumes that the differences between the two Seyfert 
types are due to geometric factors in the nuclear region. If this is so we should look for differences in the 
warm emission component between the Seyfert types. Let us therefore turn our attention to the warm dust 
component that we have suggested it is directly related to the nuclear and circumnuclear dust emission. If 
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there is a molecular torus obscuring the nuclear region and this is sufficiently thick (as the models predict, 
see e.g. Pier & Krolik 1992; Granato & Danese 1994) we shall observe differences in the warm component 
emission produced by each of the two Seyfert types. 

In order to compare the parameters obtained from the inversion of the SEDs for the different Seyfert 
types, we establish two different groups of objects within the CfA sample: Seyfert Is, including galaxies 
classified as Seyfert 1.5, and Seyfert 2s, including Seyfert 1.8 and 1.9. Therefore, the sample gets divided 
into 22 Seyfert 1 and 24 Seyfert 2. If we drop those objects not observed or those whose data do not allow 
obtaining a well defined warm component, the sample with adequate data for this analysis consists of 18 
Seyfert 1 galaxies and 22 Seyfert 2 galaxies. 

If the dust responsible for the warm emission component is located within the molecular torus that 
hides both the active nucleus and the broad line region, we expect to find differences between some of the 
parameters that define the warm component of the Seyfert 1 and 2 objects. If however the dust is located 
outside the molecular torus, then the geometry does not play anymore a key role and we should not see 
differences between the two types of Seyfert galaxies. 

6.1. Temperature differences 

Regarding the temperature distribution, it should be realized that in both Seyfert classes we are 
sampling dust emission produced under similar physical conditions, probed with the same wavelength 
range, and heated in either case by the active nucleus or by circumnuclear starforming regions. Therefore as 
temperature is an intensive quantity as opposed to extensive, we expect not to see any difference between 
the average temperature of the warm component in the two Seyfert types. 

Figure ^ shows the distribution of warm component temperature for the two groups of Seyfert galaxies. 
It is seen that the median temperature of the Seyfert Is is slightly higher than the median of the Seyfert 
2s though the differences are small, IZwl excepted. As we pointed out before, IZwl is a special object, its 
warm temperature being way far out of every other object in the sample. We have therefore dropped IZwl 
from the statistics. To understand the significance of the differences between the two distributons, we have 
applied a Kolmogorov-Smirnov (KS) test. The characteristics of the distributions are: 



(T)i =148 a = 12. 
(T} 2 =144 g = 14. 
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The result of this KS test is that there is a 78% probability that the two distributions are the same 
and so the differences between them are not significant and both groups (Seyfert 1 and Seyfert 2) enjoy as 
expected, the same temperature distribution. 

6.2. Fluxes and luminosities 

Recent studies of the mid IR emission from Seyfert galaxies claim that Seyfert 2 galaxies are weaker 
than Seyfert Is (Heckman 1995; Maiolino et al. 1995; Giuricin, Mardirossian & Mezzeti 1995; Mulchaey et 
al. 1994). This result is interpreted within the framework of the unified models as an anisotropy, resulting 
from the presence of a molecular torus with a given optical thickness in the IR. To test this claim we 
have compared the fluxes and luminosities of the warm component of the Seyfert SEDs to see if we find 
significant differences between Seyfert 1 and Seyfert 2s. 

The median values for the warm luminosities are: 

{logL warm ) x = 44.7 a = 0.3 
(logL warm ) 2 = 44.9 a = 0.8 

Statistically, the KS test shows that there is a probability 33% that both distributions are the same. 
Hence there are no significant differences between the Seyfert 1 and 2's regarding the warm emission. 

Turning now to the total far IR luminosities, i.e., the sum of the three emission components, the 
distributiones are: 

{logL IR ) 1 =45.1 (7 = 0.8 
(logL IR ) 2 = 44.9 a = 0.7 

In this case the distributions are the same with a significance level of 99% therefore it can not be 
concluded that the Seyfert Is are brighter than the Seyfert 2s in the mid and far IR. 

We now turn to search wheather there are differences between Seyfert 1 and Seyfert 2s in the ratio of 
warm component to total luminosity. First, we refer to the ratios between the warm and total fluxes. The 
mean and standard deviations of these distributions are: 



{F warm /F IR ) 1 =0.42 CT = 0.17 
(F warm /F IR ) 2 =0.28 (7 = 0.15 
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Figure |i] shows the flux ratio distributions. A KS test indicates that the distributions are different at 
the 99% level of significance. This result suggests that the Seyfert Is emit fractionally more than the Seyfert 
2s in the mid IR (warm component). This is consistent, for example, with the models of circumnuclcar tori 
of Pier & Krolik (1992;1993) and of Granato & Danese (1994), that predict anisotropy by absorption of 
nuclear emission in the mid IR. However this result should be taken with care, as it could indicate either 
that the nuclear emission is larger in the Seyfert Is relative to their total FIR flux or that the contribution 
of the host galaxy is stronger in the Seyfert 2s. To discriminate between these two possible explanations 
we must normalize the infrared fluxes with an isotropic property of the galaxies in the sample, i.e., with 
fluxes emitted at long enough wavelengths that they are not suspect of suffering extinction and thus do not 
depend on the geometry of the sources. Other authors have used [OIIIJA5000 A fluxes, hard X-ray fluxes or 
radio fluxes (see, for example, Mulchaey et al. 1994). However, the [OIII] A5000 A fluxes can be affected 
by absorption due to dust in the NLR. We prefer to use 20 cm radio emission fluxes to normalize the IR 
flux, since the radio emission is not affected by selective extinction. We have used integrated radio data 
from Edelson (1987). These data consists of VLA observations at 1.46 GHz (20 cm) with a bandwidth of 
45 MHz. The FWHM beamwidth used is 1.5 arcmin, directly comparable with our ISO data. 

The distributions of the radio normalized warm IR fluxes show the following characteristics: 



KS tests indicate that both distributions are different at a significance level of 99.9%. Therefore, the 
warm flux is indeed higher in Seyfert Is than it is in Seyfert 2s, and it can be concluded that the warm 
emission from Seyfert 2s is affected by dust extinction to a larger extent that in the Seyfert 1 galaxies. 

This result suggests that at shorter wavelengths (mid IR) the emission is still anisotropic, in agreement 
with the molecular torus models of Pier & Krolik (1992;1993) and Granato & Danese (1994). These and 
others authors have proposed different models for the absorbing material. Pier & Krolik (1992;1993), 
Granato & Danese (1994) and Efstathiou & Rowan-Robinson (1994) have modeled the absorbing structures 
as axially symmetric tori. The models proposed by Granato & Danese (1994) result in thin and extended 
tori with optical depths ranging from r w 10 to 300 in the UV band and maximum radii ranging from 
tens to hundreds of parsecs. On the contrary, the models proposed by Pier & Krolik (1993) show thin 



(log(F x 
(log(F, 



warm 



warm 



/F 2 o cm ))i =6.5 cr = 0.3 
/F 2Qcm )) 2 =6.1 a = 0.3 
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and compact accretion disks with very large optical depths with values of r ~ 1000 in the UV band, and 
compact radii with dimension of a few pc. If we consider the Seyfert 1 as canonical unobscured objects, 
and abscribe the differences found between the two Seyfert types to absorption by the obscuring torus we 
obtain a mid IR optical depth of tir « 0.4 or a tuv ~ 80 (assuming t\ oc 1/A) for the Seyfert 2 objects. 
This value is indeed very mild and within the range predicted for the thin and extended tori of Granato, 
Danese & Franceschini (1997). 

The validity of ratioing with the 20 cm radio flux has been however questioned based on the idea that 
Seyfert 2 galaxies may have more star formation in their disks than Seyfert Is (Maiolino et al. 1995). This 
would affect the radio emission, and hence the warm to radio flux ratio. It remains to test wether the 
Seyfert 2 galaxies are indeed stronger emitters of extended far IR radiation. We consider for this test the 
ratio of the cold far IR component to 20 cm radio flux. For the two groups the values are: 

(log(F cold /F 2 ocm))i =6.7 a = 0.7 
(log(F cold /F 20 c m )}2 = 6.5 a = 0.5 

These distributions are similar (60% probability) according to the KS test. Therefore, Seyfert 2s are 
not stronger emitters in extended IR radiation than Seyfert Is, and there are no reasons to suspect that 
they should be stronger radio emitters based solely on the amount of star formation occurring in Seyfert 
2s. It is also worth pointing out that the differences found between the type 1 and 2 Seyferts are restricted 
to the warm emission, while there are not differences regarding the cold and very cold emission, i.e., the 
emission from their respective galaxy disks. 

7. Summary 

We have presented far IR photometry with ISO of the CfA Seyfert sample. The data have allowed 
a detailed study of the far IR SED of these sources using a Bayesian inversion method. It has been 
shown that the mid and far IR emission of Seyfert galaxies can be explained by the emission of three 
thermal components, a warm component, associated with dust heated by the nucleus and circumnuclear 
starformation regions; a cold dust component heated by star forming region in the galaxy disk, and very 
cold dust component heated by the general interstellar radiation field. The mid to far IR output from 
Seyfert galaxies does not have a simple origin but different ingredients play an important role in it. 



- 18 - 



The comparison of cold far IR fluxes with Ha data confirms that the cold emission component that 
appeared naturally from the inversion process is related to radiation from dust in starforming regions in the 
galaxy discs. 

We find that the mid IR emission (warm component) is larger in Seyfert 1 than in 2s, suggesting the 
presence of obscuring material in Seyfert 2s. The median value obtained for the optical depth is in the 
range predicted by the thin and extended tori models. 
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Fig. 1. — Spectral Energy Distributions of the CfA Seyfert sample. Fluxes are in Janskys 

Fig. 2.— Mid and far IR SEDs of the CfA Seyfert sample. For each object, the upper panel shows the ISO 
(triangles) and IRAS (filled squares) data and the best fit to the SED obtained with the Planckian inverse 
transform. The different thermal components contributing to the fit are printed in dashed lines. The bottom 
panel shows the temperature spectrum responsible for these components 

Fig. 3. — Mid and far IR Spectral Energy Distribution of NGC 1068. As in Figure 2, ISO data are triangles 
and IRAS data are filled squares. Top: The two thermal components obtained with the inversion method 
are plotted with dashed lines. In this case, we have not considered the two first shorwavelcngth data points 
in order that the spectral range are equal to the rest of CfA sample. Bottom: Result after applying the 
inversion method to all the complete data set. 

Fig. 4. — Top: Distribution of temperatures of the warm components for the objects of the CfA sample. 
Note that the value of IZwl is clearly higher than the rest of the temperatures. Bottom: Distribution of the 
ratio between the warm and the total fluxes of the CfA sample 
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Tablc 1. ISO fluxes (in Janskys) for the CfA Seyfert sample 
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Table 1 — Continued 
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Table 2. Fluxes of NGC 1068. Fluxes are given in Janskys 
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Tabic 3. Central temperature of spectral components 
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NGC 5252 


2 


SO 


20 


56 


167 


Mrk 266 


2 


Pec. 


27 


45 


130 
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Table 3 — Continued 



Object 


Type 


Morph. type 


Ti(K) 


T 2 (K) 


T 3 (K) 


Mrk 270 


2 


SO 


18 


55 


161 


NGC 5273 


1 


SAO 


31 


63 


151 


Mrk 461 


2 


s 


20 


63 


153 


Mrk 279 


1 


SO 


24 


46 


153 


IC 4397 


2 


s? 


26 


46 


158 


NGC 5548 


1 


SAO 


23 


51 


152 


Mrk 817 


1 


s? 




45 


135 


Mrk 841 


1.5 


E 




71 


162 


NGC 5929 


2 


Sab 


28 


54 


126 


1614+35 


1.5 


S? 


19 


41 


150 


2237+07 


2 


SBa 


22 


58 


151 


NGC 7469 


1 


SABa 


27 


60 


140 


Mrk 533 


2 


SAbc 


21 


54 


149 


NGC 7682 


2 


SBab 


16 


45 


136 



Note. - 



Short wavelength data (5-11 /*m) not taken into account (See section 



4.1 



) 
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Table 4. Fluxes for each component of the SEDs (erg cm" 2 s _1 ). Fi corresponds to the very cold 
component, F 2 to the cold component and F 3 to the warm component. F tota ; is the total sum of the 

fluxes 



Object 


Type 




F 


1 






F 


2 






F 3 




Ftotal 




Mrk 334 


2 










2 


.964 


10" 


-10 


9.618 


io- 


■ 11 


3.655 


10" 


10 


Mrk 335 


1 


1. 


.326 


10" 


-11 


2 


.357 


10" 


-11 


4.211 


io- 


-11 


7.894 


10" 


11 


IZwl 


1 


1. 


135 


10" 


-10 


9 


.073 


10" 


-11 


>1.345 10~ 10 


>3.453 10- 10 


Mrk 993 


1.5 


2. 


.072 


10" 


-11 





.793 


10" 


-11 


2.322 


io- 


-11 


5.187 


10" 


11 


Mrk 573 


2 


9. 


.210 


10" 


-11 


2 


.596 


10" 


-11 


1.120 


io- 


-10 


2.300 


10" 


10 


0152+06 


2 


7. 


.836 


10" 


-11 


1 


.823 


10" 


-11 


4.062 


io- 


■ 11 


1.372 


10" 


10 


Mrk 590 2 


1 






























NGC 1068 


2 










1 


.170 


10" 


-8 


1.198 


io- 


■8 








NGC 1144 


2 


5. 


131 


10" 


-10 


1 


.702 


10" 


-10 


1.012 


io- 


-10 


7.845 


10" 


10 


Mrk 1243 


1 


2. 


.187 


10" 


-11 


1 


.209 


10" 


-11 


2.188 


io- 


-11* 


5.584 


10" 


n 


NGC 3079 


L 


2. 


.929 


10" 


-9 


1 


.344 


10" 


-9 


4.283 


io- 


10 


4.701 


10" 


9 


NGC 3227 


1.5 


5. 


.505 


10" 


-10 


3 


.850 


10" 


-10 


2.805 


io- 


-10 


1.216 


io- 


9 


NGC 3362 


2 


5. 


168 


10" 


-11 


9 


.518 


10" 


-11 


5.533 


io- 


■ 11* 


2.022 


10" 


10 


1058+45 


2 


5. 


.712 


10" 


-11 


2 


.351 


10" 


-11 


3.289 


io- 


■ 11 


1.135 


10" 


10 


NGC 3516 


1.5 










1 


.201 


10" 


-10 


1.593 


io- 


■10 


2.794 


10" 


10 


Mrk 744 


2 


1. 


923 


10" 


-10 


5 


.239 


10" 


-11 


4.684 


io- 


-li 


2.924 


10" 


10 


NGC 3982 


2 


3. 


.597 


10" 


-10 


2 


.951 


10" 


-10 


1.121 


io- 


-10 


7.669 


10" 


10 


NGC 4051 


1 


8. 


.771 


10" 


-10 


2 


.690 


10" 


-10 


3.053 


io- 


-10 


1.1.451 10~ 9 


NGC 4151 


1 


1. 


519 


10" 


-10 


2 


.961 


10" 


-10 


8.346 


io- 


■10 


1.282 


10" 


9 


NGC 4235 


1 


3. 


.132 


10" 


-12 


2 


.557 


10" 


-li 


3.950 


io- 


■ li 


6.820 


10" 


11 


Mrk 766 


1.5 


1. 


.537 


10" 


-10 


1 


.187 


10" 


-10 


1.741 


io- 


■10 


4.465 


10" 


10 


Mrk 205 


1 


3. 


669 


10" 


-li 


3 


.464 


10" 


-li 


3.344 


io- 


n 


1.047 


io- 


10 


NGC 4388 


2 


4. 


980 


10" 


-10 


4 


.286 


10" 


-10 


4.540 


io- 


10 


1.380 


10" 


9 


Mrk 231 


1 










2 


.711 


10" 


-9 


9.569 


io- 


10 


3.668 


10" 


9 


NGC 5033 


2 


1. 


.822 


10" 


-9 


4 


.121 


10" 


-10 


2.586 


io- 


10 


2.502 


10" 


9 


1335+39 


2 


8. 


.249 


10" 


-11 


2 


.323 


10" 


-11 


3.340 


io- 


■ 11 


1.391 


10" 


10 


NGC 5252 


2 


2. 


.250 


10" 


-11 


2 


.232 


10" 


-11 


2.812 


io- 


■ 11 


7.294 


10" 


11 


Mrk 266 


2 


2. 


.049 


10" 


-10 


3 


.164 


10" 


-10 


1.168 


io- 


■10 


6.381 


10" 


10 
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Tabic 4 — Continued 



Object Type Fi F 2 F 3 Ytotal 



Mrk 270 


2 


1.244 10~ n 


5.686 10~ 12 


1.565 10- 11 


3.378 10 


NGC 5273 


1 


5.174 10~ n 


3.219 10~ n 


1.576 10- 11 


9.969 10 


Mrk 461 


2 


5.384 10~ 12 


1.874 10~ n 


2.974 10~ n 


5.386 10 


Mrk 279 


1 


7.046 10~ n 


3.165 10~ n 


6.831 10~ n 


1.704 10 


IC 4397 


2 


1.194 10~ 10 


3.566 10- 11 


3.635 10~ n 


1.913 10 


NGC 5548 


1 


1.769 10- 11 


5.768 10- 11 


1.126 10~ 10 


1.898 10 


Mrk 817 


1 




1.680 10" 10 


1.357 10~ 10 


3.037 10 


Mrk 841 


1.5 




2.913 10~ n 


6,982 10- 11 


9.895 10 


NGC 5929 


2 


5.487 10" 10 


2.282 10~ 10 


1.654 10~ 10 


9.423 10 


1614+35 


1.5 


3.392 10" 11 


2.227 10~ n 


1.923 10~ n 


7.542 10 


2237+07 


2 


1.297 10~ n 


3.840 10~ n 


5.347 10~ n 


1.048 10 


NGC 7469 


1 


1.231 10~ 9 


9.128 10~ 10 


4.390 10~ 10 


2.583 10 


Mrk 533 


2 


1.532 10~ 10 


2.594 10~ 10 


2.217 10~ 10 


6.643 10 


NGC 7682 


2 


1.058 10- 11 


1.796 10- 11 


4.352 10~ n 


7.206 10 
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Table 5. Luminosities for each component of the SEDs (10 44 erg s Li corresponds to the very cold 
component, L 2 to the cold component and L 3 to the warm component. Ljr is the total luminosity 



Obiect 


Tvdc 


- Lj l 






*- i IR 


Mrk 334 


2 




25.19 


9.00 


34.19 


Mrk 335 


1 


1.72 


3.06 


5.47 


10.25 


IZwl 


1 


81.49 


65.14 


>96.57 


>243.20 


Mrk 993 


1.5 


0.94 


0.36 


1.06 


3.76 


Mrk 573 


2 


5.31 


1.50 


6.46 


13.27 


0152+06 


2 


4.57 


1.06 


2.37 


8.10 


Mrk 590 2 


1 










NGC 1068 


2 




30.06 


31.93 


61.99 


NGC 1144 


2 


82.51 


27.37 


16.27 


126.15 


Mrk 1243 


1 


5.30 


2.93 


5.30 


13.53 


NGC 3079 


L 


7.68 


3.52 


1.12 


12.32 


NGC 3227 


1.5 


1.52 


1.06 


0.77 


3.35 


NGC 3362 


2 


5.15 


9.48 


5.51 


20.14 


1058+45 


2 


9.38 


3.86 


5.40 


18.64 


NGC 3516 


1.5 




1.67 


2.21 


3.88 


Mrk 744 


2 


3.06 


0.83 


0.75 


4.64 


NGC 3982 


2 


0.62 


0.51 


0.19 


1.32 


NGC 4051 


1 


0.81 


0.25 


0.28 


1.34 


ATP P A 1 ^ 1 

lNvj\_y 4101 


1 
1 


U.ZD 


U.Ol 


1 AA 
1 .44 


9 91 
Z.Z1 


NGC 4235 


1 


0.035 


0.29 


0.45 


0.77 


Mrk 766 


1.5 


4.84 


3.74 


5.49 


14.07 


Mrk 205 


1 


35.53 


33.55 


32.39 


101.47 


NGC 4388 


2 


6.12 


5.26 


5.58 


16.96 


Mrk 231 


1 




888.77 


313.71 


1202.48 


NGC 5033 


2 


3.14 


0.71 


0.44 


4.29 


1335+39 


2 


6.43 


18.12 


26.05 


50.60 


NGC 5252 


2 


23.20 


23.02 


29.01 


75.23 


Mrk 266 


2 


30.02 


46.37 


17.11 


93.50 
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Tabic 5 — Continued 



Object 


Type 


Li 


L2 


L3 


Ur 


Mrk 270 


2 


0.19 


0.09 


0.24 


0.52 


NGC 5273 


1 


0.13 


0.08 


0.04 


0.25 


Mrk 461 


2 


0.27 


0.92 


1.47 


2.66 


Mrk 279 


1 


12.63 


5.68 


12.25 


30.56 


IC 4397 


2 


4.97 


1.48 


1.51 


7.96 


NGC 5548 


1 


0.94 


3.06 


5.98 


9.98 


Mrk 817 


1 




32.15 


25.97 


58.12 


Mrk 841 


1.5 




7.51 


18.00 


25.51 


NGC 5929 


2 


7.25 


3.02 


2.19 


12.46 


1614+35 


1.5 


5.15 


3.38 


2.92 


11.45 


2237+07 


2 


1.57 


4.64 


6.47 


12.68 


NGC 7469 


1 


60.72 


45.02 


21.65 


127.39 


Mrk 533 


2 


24.81 


42.01 


35.90 


102.72 


NGC 7682 


2 


0.59 


1.00 


2.42 


4.01 



Mrk 334 



Mrk 335 



100.00 



3) 



1.00 



0.10 



100.00 



3> 



1.00 



0.10 



0.01 



\( lira) 
Mrk 573 



10.00 r 



0.01 




10 



100 



Mrk 993 



' t 1 



0.01 



10 



100 



\( fira) 
0152+ 06 




t t 

t < H 



t t 



0.01 .... I 
10 



100 



~h(fim) 



\(fj,m) 



Mrk 590 



NGC 1144 



0.01 



10 



100 



Mrk 1243 



\( fj.rn) 
NGC 3227 



100.00 



3> 



1.00 



0.10 





t M t+ 



0.01 



10 



100 



NGC 3079 



t t 



0.01 



10 



100 



\( firn) 
NGC 3362 



ft +++ 



10 



100 



\(jj,m) 



1058+ 45 



NGC 3516 



ft \ 



0.01 



10 



100 



Mrk 744 



\( lira) 
NGC 4051 



100.00 



3> 



1.00 



0.10 





t t 



0.01 



10 



100 



NGC 3982 



0.01 



10 



100 



\( fira) 
NGC 4151 



10 



100 



\(jj,m) 



NGC 4235 



Mrk 766 



t t 



0.01 



10 



100 



Mrk 205 



\( fj.rn) 
Mrk 231 



100.00 





t t 



0.10 



0.01 



10 



100 



NGC 4388 



0.01 



10 



100 



\( firn) 
NGC 5033 



t t 



tt 



10 



100 



\(jj,m) 



1335+ 39 



NGC 5252 



0.01 



10 



100 



Mrk 266 



Xf firn) 
NGC 5273 



100.00 



3> 



1.00 



0.10 





0.0? 



10 



100 



Mrk 270 



0.01 



10 



100 



Mrk 461 



< t III 



/o 



7 00 



\(jj,m) 



1614+ 35 



2237+ 07 



0.01 



10 



100 



NGC 7469 



\( firn) 
NGC 7682 



100.00 



3> 



1.00 



0.10 





0.01 



10 



100 



Mrk 533 



+ t + + 



0.01 



10 



100 
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NGC 1068 
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